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One of the confident predictions of QCD is that at sufficiently high temperature and/or density, 
hadronic matter should undergo a thermodynamic phase transition to a colour deconfined state 
of matter - popularly called the Quark-Gluon Plasma (QGP). In low energy effective theories of 
Quantum Chromodynamics (QCD), one usually talks of the chiral transition for which a well defined 
order parameter exists. We investigate the dissociation of pions and kaons in a medium of hot quark 
matter decsribed by the Nambu - Jona Lasinio (NJL) model. The decay widths of pion and kaon 
are found to be large but finite at temperature much higher than the critical temperature for the 
chiral (or deconfinement) transition, the kaon decay width being much larger. Thus pions and even 
kaons (with a lower density compared to pions) may coexist with quarks and gluons at such high 
temperatures. On the basis of such premises, we investigate the process of hadronization in quark- 
gluon plasma with special emphasis on whether such processes shed any light on acceptable low 
energy effective theories of QCD. 
Quantum Chromodynamics (QCD) is believed to be the underlying theory of strong interactions. Although enor- 

• mously successful, inasmuch as it is part of the standard model of physical interactions, the applicability of QCD is 
studying nuclear interactions or hadronic processes at low energies is still limited. The difficulties are largely technical, 
associated with the non-perturbative features of the theory and hence, a large amount of effort is being devoted to 

jy-j ' finding effective lagrangians which contain features compatible with the low energy limit of QCD. This is the thrust of 
the present workshop. The philosophy of my talk is somewhat complementary; I am going to discuss QCD phenomena 
which occur in very energetic nuclear collisions, where one hopes perturbative QCD does play an important role, at 
, least during the initial time period. 

A confident prediction of QCD is that at very high temperature and/or density, the bulk properties of strongly 
interacting matter would be governed by the coloured QCD degrees of freedom - quarks and gluons- rather than the 
usual hadrons. Such a phase is called quark gluon plasma (QGP) (!]] in the literature. Conditions conducive to the 
' formation of QGP may have existed in the early universe during the first few microseconds after the Big Bang. Also, 
such conditions may be transiently created in highly energetic collsions of large nuclei and the search for such a novel 
O ' phase of matter constitutes a major area of current research in the field of high energy physics. 

Whether QGP is separated from the hadronic world by an actual thermodynamic phase transition is an open 

• question. It has been widely postulated that such a phase transition may indeed occur, where the quarks and gluons 
convert into colourless hadrons. Recent results, showing the lack of thermodynamic equilibrium [|| in the quark-gluon 
phase in ultrarclctivistic heavy ion collisions, indicate however that such an ideal situation is unlikely. It should also 
be noted that although the persistence of non-perturbative effects till very high temperatures was suggested in the 
literature quite early on [||, it is only recently that the lattice results have confirmed that non-perturbative hadron 
like excitations could survive at temperatures far above the chiral phase transition temperature [Q . The lattice result 
for pion screening mass has been studied in ref. || . The analysis of ^ has been contradicted by Boyd et al. in ref. [|| . 
The conclusion of these authors Q is consistent with the existence of free quarks at high temperatures. On the other 
hand, Shuryak argued in a subsequent work that the non-perturbative modes, especially pion- like excitations, 
could indeed survive till temperatures above T c . Furthermore, similar results for pion screening masses are obtained 
in a- model as well It is thus imperative to understand the behaviour of such hadronic resonances, their formation, 
stability and so on, in a quark gluon medium at high temperature. We confine our attention to the case of pions and 
kaons only; these, being lighter than other hadrons, account for the bulk of the multiplicity. 

Formation of light mesons like pions and kaons, a bound state of light relativistic quarks, is an extremely difficult 
problem to handle in QCD, where all the troublesome features of non-perturbative QCD appear. We therefore employ 
the usual practice of looking at the pion and kaon as Goldstone bosons arising from the spontaneous breaking of the 
chiral symmetry, a feature most suitably addressed in the Nambu Jona-Lasinio (NJL) model Q . 

The formulation of NJL model in flavour SU(3) was first introduced by Hatsuda et al. || and Bernard et al. [fl0|| . 
The three flavour NJL model Lagrangian is written in terms of u, d and s quarks, the interaction between them being 
constrained by the SU(3)l <8> SU(3)r chiral symmetry, explicit symmetry breaking due to the current quark masses 
and the U(1)a breaking due to the axial anamoly H. The full Lagrangian with KMT (Kobayashi- Maskawa -'t-Hooft) 
term is given below b| . 
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+ g D [detq~i(l - 75)^ + h.c] (1) 

where the quark fields qi has three colours (N c = 3) and three flavours (Nt = 3) and A a (a = 1, 8) are the Gell-Mann 
matrices. The quark mass matrix is given by m = diag(m Ul m c i, m s ). 

In the mean field approximation, the quark condensates at finite temperature are given by, 

« Mi »= ~2N C J 0^/(E( P ) (2) 

where Ei p is the quark single particle energy for the i-th specie and f(Ei P ) = 1 — rii p — hi pi rii p and Ui p being 
the Fermi-Dirac distributions for quarks and anti-quarks, respectively. If quark chemical potential is zero, then 



= [exp(E ip /T) + 1] 
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The temperature dependent constituent quark masses Mi are obtained from the expressions below, 

(3) 

where 

«M»Ea, « dd »=/?, << ss >>= 7 (4) 

We now want to investigate the decay of pionic and kaonic excitations, the properties of which we assume to be 
given by the NJL model. It should be mentioned here that at temperatures above the critical temperature, these 
mesonic excitations are more like resonances with large effective masses In the following, we study the 

decay width of such pseudoscalar excitations in the hot quark medium as a function of temperature, starting with 
the Lagrangian given above in equation (Fjj) . 

The quark mass Mi appearing in eq. ||) and in eq. (||) is a very important ingredient in our calculation. In the 
absence of any medium and/or dynamic effect, Mi is the current quark mass. On the other hand, we know that due 
to the spontaneous breakdown of the chiral symmetry, quarks attain the value of the constituent quark mass fl2|| . 

In the present calculation we have used the parametrisation of rcf. jlO| ( A = 631.4, <? S A 2 = 3.67, gf£>A 5 = -9.29 
and current mass m Ui <i(rn s ) = 5.5 (135.7) MeV ) to calculate the quark and meson masses. The constituent quark 
masses are calculated using the gap equations (eq. ||). These quark masses are then put into the dispersion equation 
for mesons to get dynamical masses of mesons (n and K, here). 

1 + 2G^\ij{oj,q^G) = (5) 

where Ily is the one loop polarization due to u and d quark for pions and u or d and s quark for kaons. is the 
coupling constant with <f> coerresponding to 7r or K. The general expression for polarization function is 
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where N c is the number of colours and k = p + q. The energies E p = Vp 2 + Mi 2 and E k = yj (p + q) 2 + M 2 2 . For 
pion, Mi = M 2 = M u . For kaon, Mi = M u u\ and M2 — M s . n k and n p are the Fermi-Dirac distribution functions 
defined earlier. The pseudoscalar couplings are, 

G-n = g s + g D l 
G K ± = 9s + 9d(3 

G K ° = 9s + 9dol (7) 
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where a, (3 and 7 are denned in eq. (ji]). 

The decay width is evaluated using the imaginary part of the eq. 



as given below, 



I\4 = 



(8) 



where G<t, q is the empirical meson-quark coupling as obtained in NJL. Here we have used G vq = 3.5 and Gxq = 3.6 
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Figure 1 : Variation of quark and meson masses with temperature. 
Figure 2 : Variation of meson decay widths with temperature. 



The variation of quark and meson masses is shown in figure 1. The u or d quark masses starting from 135 MeV 
drops to the current quark mass value just after a temperature of 200 MeV. On the other hand, the drop in the 
strange quark mass is much smaller around that temperature, showing the effect of explicitly broken chiral symmetry, 
by a larger amount, in the SU(3) sector. Pion and kaon both show a similar qualitative behaviour. The masses of 
pion and kaon remain constant at their free masses upto a temperature of 200 MeV but increases sharply after that, 
the pion mass rising to 900 MeV and kaon mass to 1000 MeV around a temperature of 450 MeV, thus giving a slower 
increment for kaons compared to pions.The difference in the behaviour of kaon and pion can be attributed to the 
difference in u and s quark masses. 

Figure 2 shows that the decay width, for both pion and kaon, is very high at high temperature and decreases 
with decreasing temperature, going to zero at around T — 0.2 GeV. It is worth noticing that at around the same 
temperature, the effective pion mass attains the value of the free pion mass (figure 1). The decay width of kaon is 
around 3 GeV where as that of pion is around 1.4 GeV at T = 500 MeV. This is very significant for two reasons. 
Firstly, our results show that though there will be pions and kaons along with the quarks at high temperature phase, 
the numbers of mesons will be very small due to their large decay width. Moreover, the number of kaons will be 
much less compared to pions at high temperature phase, though both the mesons will become stable around the same 
temperature (below 200 MeV). 

We thus find that even without any consideration of the detailed evolution and dynamics of the system, the mesonic 
modes in a hot quark medium are found to become important around a temperature of 200 MeV. Though the question 
whether this is a signature of a phase transition cannot be addressed within the present framework, the fact that most 
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of the pions and kaons decay into quarks, owing to a large decay width at temperatures higher than T = 200 MeV, is 
a remarkable finding. Moreover, both the pionic as well as kaonic modes start becoming important at about the same 
temperature, thus providing a hint of some kind of a transition temperature. It is therefore tempting to attempt a 
microscopic investigation of the process of hadronization within the NJL model. 

Recently there have been some attempts to study the formation of hadrons in quark matter using different semi- 
microscopic approaches. These studies can be characterized either as model dependent calculations or the 
computer codes based on the string phenomenology , or other phenomenological description of hadronization [[L5j . 
None of these approaches account for the essential lack of equilibrium in the quark-gluon phase. Some efforts have 
also been made to estimate hadronization within the parton cascade model by introducing a cut-off to mimic the 
non-perturbative effects To the best of our knowledge, the first study aimed at investigating the dynamical 
process of hadronization in a non-equilibrated quark-gluon system from a physically transparent approach was in [ |L7| . 

We have earlier shown || that perturbative estimates of the gluon-gluon, quark-gluon and quark-quark cross sections 
allow us to study the evolution of the quark-gluon matter formed in ultrarelativistic heavy ion collisions by visualizing 
the quarks as Brownian particles in a hot gluonic thermal bath. Let us start with such a premise, which, in our 
opinion, describes the non-equilibrium aspects of the evolution in a physically transparent manner. 

Hadronization in such a system can be studied in the light of Smoluchowski's theory of coagulation in colloids, 
which was elaborated further by Chandrasekhar [H . This theory suggests that coagulation results as a consequence 
of each colloidal particle being surrounded by a sphere of influence of a certain radius R such that the Brownian 
motion of a particle proceeds unaffected only so long as no other particle comes within its sphere of influence. When 
another Brownian particle does come within a distance R of the test particle, they form a two body cluster. This 
cluster also describes a Brownian motion but at a reduced rate due to its increased size/mass. The process continues 
further till a single cluster of all the particles is formed. 

In order for this stochastic scenario of cluster formation to apply to a system of Brownian quarks in the hot gluon 
bath, it is essential that each quark has an appropriate sphere of influence of radius r. Obviously, this radius r will 
depend on the spin-isospin combination of final cluster (whether the final cluster is scalar, pseudoscalar, vector or 
axial vector meson or even a baryon). Mesons are formed when one quark and one antiquark with proper quantum 
numbers come within the spheres of influence of each other. (Clusters with greater numbers of quarks and antiquarks 
(e.g. baryons) can also be formed by imposing the conditions of colour neutrality and charge balance properly.) 
This implies that the radius of the sphere of influence corresponds to the correlation length between the quarks in 
the proper hadronic channel. In other words, this is the screening length of the corresponding hadrons in the hot 
quark-gluon matter. (One can immediately see, without further ado, that the rate of pion formation in the hot quark 
matter should be rather small at high temperatures and increase with falling temperature as the pion screening length 
(inverse of the screening mass) decreases with increasing temperature in all realistic pictures.) The pions formed at 
very high temperatures are most likely to decay back into quarks and antiquarks, because of the large decay width at 
high temperatures. We adopt, for the sake of consistency, the NJL model estimates for the pion screening mass Jl^] 
in the present work. 

The rate of pion formation from Brownian quarks as a stochastic process, as is evident from the preceding discussion, 
depends on the number of quarks (antiquarks) falling into the sphere of influence of another antiquark (quark), the 
number of pions decaying back to quarks and antiquarks and also the change in the pion density due to the expansion 
of the system. Thus the rate equations can be written as, 

''" r " n m n q A^ <v-f> -T™ al - ^ (9) 



^- = ^{n u nu + n d n 3 )A7Tr 2 <v-f> -T tota \ ^ ul - l(d<i) - ^ (10) 



dn rr« 1,_ _ , _ _ N „_ r 2 < > _ T total "tt° 



_|i = T total g ^ + T total gg ^ qt - n qt n q Anr 2 <v-f> +T total w ^ j - (11) 

In eqs.(pPlO|), the first term is the rate of pion formation (a = + or -); the second term is the rate of pions decaying 
back to quarks and the third term is due to Bjorken (longitudinal) expansion of the system. stands for u or d 
(we ignore s and other heavier flavours). < v ■ f > (the average relative velocity in the radial direction) is calculated 
using the Jiittner distribution, 

f(x,p) = e-zW*) (12) 



There would also be a corresponding rate equation for the antiquarks, which looks exactly like eq. ( |11| ) and hence 
not explicitly written. In eq. ( |ll"| ) the T total g ^ qq as well as Y total gg -> qq stand for the corresponding net quantities. 

As mentioned earlier, we are considering a non-equilibrated quark matter and hence the pions formed will also 
be out of equilibrium. This is taken into account by multiplying the relevant distribution functions with the ratios 
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r q = n q /n eq and rv = n^/n^ where n q and n eq are non-equilibrium and equilibrium densities of quarks and n„ and 
n e7r are non-equilibrium and equilibrium densities for pions. 

In all these expressions, the appropriate masses are the effective masses including the current as well as thermal 
contribution, , whose importance in determining the dynamics of the hot quark matter has been well established. For 
quarks (antiquarks), this is 



l eff 



\J m q (curr) 2 + m q (thermaiy 



where fB0U2 



m t (thermal) = (1 



2 A 3 ; 



and m q (curr) is taken to be 10 MeV. For gluons the thermal mass is, 

m g (thermal) — —g s T 



(13) 



(14) 



The running coupling constant a s as a function of temperature is given by 

12tt 



(33 - 2n f )ln 



A 2 



(15) 



with Q 2 = m eff 2 (T) + 9T 2 . 

Simultaneously, we must take account of energy momentum conservation which, for a Bjorken flow, corresponds to 
the following equation 



de 
di 



P 



(16) 



where e = e to tai = e g + 1 
velocity of sound, as in 



We also include the one loop correction to e„ j23|. e and P are related through the 
For a complete description of the system, eqs. (g), (f[o|), ( |ll| ) and jT6| ) must be solved 



self-consistently. The initial conditions are taken from |2j for RHIC energies. The initial time (t g ) is the time when 
gluons thermalise (=0.3 fm), where r q = 0.15, r g =l and r w is taken to be 0. The temperature at this time is 500 
MeV. Note that we are working at y = so that t and r are the same and the baryon chemical potential is zero. 
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Figure 3 : Time evolution of quark and pion number densities for various values of QCD parameter A. 
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Figure 3 shows the variation of pion and quark number densities with time, for various values of the QCD parameter 
A. In all three cases, we find the same qualitative feature that pions start appearing in the system quite early on but 
they become appreciable in number only after some time. At late times the system is dominated by pions. This cross 
over occurs at t > 4 fm for A = 0.2 or 0.3 GeV while for A= 0.4 GeV this happens at t ~ 6 fm. 

Figure 4 shows the variation of temperature with time. Obviously, there is a dramatic effect of the QCD parameter 
A. In all the cases, there is a change at T ~ 215 MeV, corresponding to t ~ 3.5 fm; the variation of temperature with 
time becomes slower, as is expected in the mixed phase of a first order phase transition. At A=0.2 GeV, this occurs 
for a very short period of time, before the system starts cooling again. The duration of the constant temperature 
configuration increases with A, and for A=0.4 GeV, it persist upto 9 fm before the temperature of the system starts 
falling again. 
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Figure 4 : Time evolution of the Temperature for various values of QCD parameter A. 

Obviously, this is a clear indication of an apparent first order transition. Microscopically, the appearance of the 
mixed phase at a temperature of ~ 215 MeV can be understood from the fact that the pion decay width goes to 
zero at such a temperature. All the pions that were formed earlier in the system tended to decay back to quarks 
and antiquarks on a fast time scale. Only after the pion decay width becomes small would the formed pions become 
stable. 

The apparently desirable features of the dynamical deconfining transition seem to arise quite naturally in this 
scenario. It is therefore important to test if all the conservation laws are obeyed during the evolution. To this end, 
let us check if entropy increases steadily during the entire evolution. 
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Figure 5 : Time evolution of entropy with time (A = 0.3 GeV). 
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The decrease in entropy occurs precisely where the gluons drop out of the system. This, in hindsight, was to be 
expected, since the sudden decrease in the number of degrees of freedom in the system should lead to a decrease in 
entropy. The NJL model, in the present form, does not account for the gluonic degrees of freedom. It thus appears to 
us that in order to obtain an effective low energy lagrangian from QCD capable of describing the dynamical evolution 
of the quark-gluon system, the role of the gluons must be taken into account. Incorporation of the dilaton field to 
account for the QCD trace anomaly in the NJL model |24|] is indeed a promising step in this respect, but the entropy 
carried by the gluons at high temperature should also be accommodated. Such a study is on our current agenda. 

It is a great pleasure to thank the organisers of the workshop, and especially Prof. J. da Providencia, for the 
kind invitation and their warm hospitality. This talk is based on a continuing collaboration the author has had with 
Jan-e Alam, Abhijit Bhattacharyya, Sanjay Kumar Ghosh and Bikash Sinha over the past several years. I take this 
opportunity to thank them all. 
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